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Abstract 

We present Hartree-Fock-Bogoliubov (HFB) calculations of the 
ground states of even Mg isotopes. A Skyrme force is used in the mean 
field channel and a density-dependent zero-range force in the pairing 
channel. 40 Mg and 20 Mg are predicted to be at the two-neutron and 
two proton drip-lines respectively. A detailed study of the quadrupole 
deformation properties of all the isotopes shows that the ground states 
of 36 > 38 ' 40 Mg are strongly deformed with significantly different defor- 
mations for the neutrons and protons. Our study supports the disap- 
pearance of the N = 28 shell gap in the Mg and Si isotopes. 
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1 Introduction 

The existence of deformed drip-line nuclei has been recently questionned [1] 
on the basis of a detailed analysis of the available data. Tanihata et al. have 
argued that with few exceptions, the last bound isotopes of a chain could be 
related to a sub-shell closure in a spherical shell model. The number of cases 
supporting their analysis is however limited since the neutron drip line is 
known only up to Z=9. A justification relying on results of fully microscopic 
calculations is not easy to obtain. For evident reasons of simplicity, most 
studies of drip-line nuclei have been restricted to spherical shapes [2, 3, 4, 
5, 6]. /,From their outcomes, it has been suggested that the single particle 
gaps associated with magic numbers significantly decrease near the drip lines 
[2]. Recently, two recent studies of deformation in the Mg[7] and Si-Ca[8] 
regions in which pairing correlations were described within the constant gap 
BCS approximation have found sizable deformations for some drip line nuclei. 
However, no final conclusions can be drawn from these studies since the BCS 
approximation is known to break down when bound and unbound single- 
particle states interact via pairing correlations [9]. 

In a recent publication [10], we have presented a method based on the 
solution of the Hartree-Fock-Bogoliubov (HFB) equations on a 3-dimensional 
cartesian mesh. It describes consistently both the pairing correlations and 
the asymptotic behaviour of the wave functions of weakly bound systems. 
In this calculation, the effective interaction is considered separately in two 
different channels: a Skyrme force to describe the mean-field (particle-hole 
channel) and a density-dependent zero-range force for the pairing correlations 
(particle-particle channel). Tests of the method on the even Ni isotopes have 
shown that it accounts correctly for the interactions between bound and 
continuum single-particle states. Since all the Ni isotopes are spherical, the 
analysis of the results were easier because of the several degeneracies of single 
particle and quasi-particle (qp) energies. On the other hand, the potential 
of investigation of the method is broader than shown by the application 
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presented in ref.[10]. Indeed, it can as well handle deformed nuclei since all 
single particle states are calculated within a three-dimensional geometry and 
allows the description of triaxilly deformed nuclei. 

This work is mostly devoted to a study of Mg isotopes. This isotopic 
chain has been selected because large deformations have been predicted for 
the heaviest isotopes on the basis of the BCS approximation. Also, together 
with neighboring isotopic chains, these nuclei have been extensively studied 
experimentally. Many experimental works have looked for deformation effects 
at iV = 20 in the A = 30 mass region. Results on isotope shifts of Na 
isotopes[ll, 12] are consistent with an appearance of deformations at N = 20 
although not yet fully conclusive because the measure of the quadrupole 
moment of 31 Na is still out of reach experimentally. More convincing evidence 
comes from measurements of masses in 31 ' 32 Mg [13], 30_32 Na and 31 ~ 33 Mg[14], 
although is not exempt of ambiguities[15]. However, the flattening of the 
curve in Mg isotopes around N = 20 appears now to be well established [16]. 
Systematics of the energies of the first 2 + states [17] may also be interpreted 
as a sign of deformation in 32 Mg. Finally, a large B(E2) value[18] has recently 
been measured for the transition between the first 2 + state and the ground 
state in 32 Mg. Most experimental data point towards the existence of a 
spherical shell closure at N = 20 for the heavier elements: 33 Al [19] and 
33-35 Si [20, 21, 22]. On the other hand, the vicinity of the one-proton drip 
line has been explored by Langevin et al. [23] by the observation 23 Si, 27 S, 
31 Ar and 35 Ca. Recently Suzuki et al. [24] have deduced matter rms radii 
for a series of Na and Mg isotopes from interaction cross sections. 

Several shell-model calculations [25, 26, 27, 28, 29] have been performed 
for Mg isotopes. They concluded that the inclusion of the f shell is crucial 
for a correct reproduction of the properties of 32 Mg. These studies also 
support the occurence of prolate quadrupole deformation in 32 Mg. Until 
very recently[30], systematic mean-field calculations in Na and Mg isotopes 
have not been numerous. The present work contains a drip-line-to-drip- 
line investigation of the deformation properties of the even Mg isotopes. In 
addition, at the neutron drip line, we have also studied the neighbouring 
even Z elements (Ne and Si). 
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2 Strength of the pairing interaction 



In our method, distinct parametrizations for the interaction are chosen for the 
mean-field and the pairing channels. For the particle-hole channel, a Skyrme 
force is used. We have studied two well established parameter sets. The 
first one is SIII [31], which has been extensively and often successfully used 
to study deformation properties. The second one, Sly4 [32] is more recent. 
It has been constructed with a special attention paid to the properties of 
nuclear matter and neutron matter. It is therefore expected to provide a 
more realistic isospin dependence than previous forces. With such a choice 
of two parametrizations of the effective force which are known to work equally 
well in the region of stable nuclei should provide a measure of the uncertainity 
of microscopic effective theory at large T z . 

For the pairing channel, a density-dependent zero-range force is used: 

\/ P (r 1 ,r 2 )=\/ (l-P (T )(l-^)5(r 1 -r 2 ) , (1) 

Pc 

where P a is the spin exchange operator and p(r) the total nuclear density. V 
is the strength of the force, chosen to be the same for neutrons and protons, 
and p c is a constant which determines the density dependence. In this work 
we take p c equal to 0.16fm~ 3 . With such a value close to the nuclear satu- 
ration density, the pairing force is strongly attractive at the nuclear surface. 
This property has been discussed in [4] (see also ref[33]). A smooth cut-off 
of the pairing interaction at an energy of 5 MeV above the fermi level has 
been introduced following the procedure explained in ref[34]. 

In ref. [10], we have shown that our numerical method provides reliable 
results for the properties of drip line nuclei when two numerial conditions 
are satisfied. Firstly, the dimension of the box in which the HF equations 
are solved must be large enough to describe correctly the tail of the wave 
functions. Secondly, since the imaginary time method used to solve the 
mean-field equations can only be applied to a limited number of orbitals, the 
number of these orbitals must be large enough to ensure a correct description 
of the continuum wave functions ; especially those associated with the low- 
lying resonance states. In the present application to Mg isotopes, we have 
used values for the box size (15 fm) and for the number of single-particle 
wave functions (70 for neutrons and 35 for protons) that have been shown 
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to be sufficient in the calculation of Ni isotopes. Since this study is devoted 
to much lighter nuclei, we are confident that these values are adequate. We 
have also checked this point with tests on a selection of nuclei, which have 
confirmed the conclusions of the previous study. 

The Mg isotopes are not an optimal region of the nuclear chart for an 
adjustment of the pairing strength Vq. Indeed, due to the weakness of pair- 
ing correlations in light nuclei, the determination of experimental qp energies 
from binding energies is not very accurate. Moreover, the variation of these 
binding energies is affected by the strong dependence of the shape of Mg 
isotopes on the neutron number, which must be taken into account in the 
process determining the qp energies. Previous studies of superdeformation 
in both the A=150 and 190 mass regions have shown that the value Vq = 
lOOOMeVfm 3 for SIII and SLy4 gives results of good quality. We have there- 
fore investigated whether one could use these values for the Mg isotopes. To 
this end, we have compared the experimental and theoretical values obtained 
for the two-neutron separation energies S 2n (N, Z), calculated from the dif- 
ference between the HFB ground state energies of neighbouring even nuclei. 
The quantity — 2A n , where A n is the HFB chemical potential, provides another 
approximation of S 2n (N, Z). In order to investigate the sensitivity of results 
to V , we have repeated some calculations with the value V r =700MeVfm 3 . 
Results are plotted on Fig. 1 together with the experimental data. The S 2n 
values are almost insensitive to Vq. They are also rather similar for both 
Skyrme parametrizations. The general features of the experimental data are 
satisfactorily reproduced. Experimentally, the shell effect at N=20 disap- 
pears in 32 Mg in agreement with the large scale shell model calculations of 
ref [29]. In this respect, the results obtained with SIII are more consistent 
with the experimental data than those obtained with SLy4, although in both 
cases the shell effect is weak. The problems associated with the description 
of 32 Mg will be discussed in a forthcoming publication. 

since the strength V = lOOOMeVfm 3 gives a reasonnable agreement with 
the experimental data for the S 2jl values, we have decided to use it for both 
SIII and SLy4. We have also checked that it leads to qp energies in 24 Mg 
close to the experimental ones. More importantly, the fact that the same 
value can be used for light nuclei and for heavy superdeformed ones is a 
first indication that a unique parametrization of a density-dependent zero- 



5 



range pairing interaction may be valid in the whole nuclear chart. With the 
adopted strength, pairing correlations vanish for protons in 36_40 Mg and for 
neutrons in 20 Mg and 26 Mg. 

On Fig. 2, we compare the experimental data for S 2jl as a function of A 
[16] to the prediction of different models. In our calculations the two-neutron 
drip line is located between 40 Mg and 42 Mg for both SIII and SLy4. The 
lightest Mg isotope predicted to be bound against two proton decay is 20 Mg. 
Its two-proton separation energy calculated with SIII is equal to 3.31 MeV; 
a number that is to be compared with the experimental value 2.33 MeV. 
Other theoretical approaches designed to reproduce nuclear masses lead to 
a better agreement with data than our calculation. However, outside the 
experimentally known region, their predictions for S 2n show irregularities 
the interpretation of which is difficult. Nevertheless, the trends obtained 
with the relativistic mean field (RMF) theory [7] and the finite-range droplet 
model (FRDM) [35], are similar except for 20,22 Mg and for a slight increase of 
S 2n from 38 Mg to 40 Mg in the FRDM. Thus it seems that the location of the 
two-neutron drip line does not depend very much on the theoretical approach, 
at least for the three types of methods we have just quoted. Hereafter we 
shall discuss the results obtained with SIII unless otherwise specified. 

3 Deformation properties of the Mg isotopes 

Fig. 3 shows the variation of the energy of the stable Mg isotopes as a function 
of their axial quadrupole moment. As expected, the energy curve of 20 Mg 
displays the well defined spherical minimum that one expects for the magic 
neutron number N=8. In light isotopes, there is a competition between oblate 
and prolate deformations. All deformed isotopes turn out prolate with the 
exception of 26 Mg. In the heavier isotopes, a deformed shell effect appears. 
It is already visible on the deformation energy curve of 32 Mg which displays 
an inflection point at Q—l. 5b. We find well defined prolate minima from 
36 Mg to the drip line nucleus 40 Mg. In these three cases, a very shallow 
oblate minimum appears at an excitation energy of 1.0 MeV. We therefore 
find that the neutron shell effect at N = 28 is suppressed at the drip line. 
/From our calculations, it turns out that the onset of deformations at iV=24 
only generates a drop of about lMeV in the curve of the evolution of S 2n 
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versus the neutron number. 

Since 40 Mg is bound by less than 2.0 MeV, we have checked if the large 
quadrupole moment of the ground state was not caused by a neutron halo 
which our calculation may not describe properly. To test the numerical 
quality of our results concerning the space distribution of outer neutrons, we 
have verified whether the number of neutrons outside a sphere of large radius 
(r >15 fm) 

iVout = / Pn(r) dh , (2) 

was stable against the box size R (p n (r) is the neutron density). We find 
N ont to be 0.430 xl(T 2 and 0.683 xl0~ 2 for box sizes R = 16 and 18 fm, 
respectively. These small and almost constant values indicate the absence 
of a halo in 40 Mg. In order to determine which part of the nuclear density 
is responsible for the calculated deformation, we have introduced quantities 
Pr( r m) (t = n,p) depending on a cutoff variable r c : 

(3 T (r c ) = (n/5)^S^ . (3) 
r z T {r c ) 

The quantities r 2 (r c ) and Q T (r c ) defined by 

^(r c ) = / rV(r) d 3 r , Q~ T {r c ) = f Q{v)p T {r) d 3 r , (4) 

Jr<r c Jr<r c 

are respectively the square radius and the quadrupole moment values of the 
fraction of the nuclear density inside the sphere of radius r c . ^From the 
results shown on Fig. 4 for 40 Mg one sees that the asymptotic values of 
the deformations (5 T are reached when r in is equal to 6.0 fm. This value is 
precisely equal to the 5/3rd of the root-mean-square radius of the nucleus 
(3.6 fm) and corresponds therefore to the sharp edge liquid drop value of the 
outer radius. This demonstrates that the neutron and proton deformation, 
j3 n and f3 p , reflect a property of the distributions of core nucleons. 

The neutron and proton parameters f3 n and (3 P and the total quadrupole 
moments are plotted in Fig. 5 for all the Mg isotopes. Values of (3 and of 
the quadrupole moments of secondary minima are also indicated. Due to 
the conjunction of the N = Z =12 deformed shell effects, the nucleus 24 Mg 
is the most deformed of the isotope chain. In our calculation, the spherical 
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shell closure at N=20 is strong enough to compensate the effect of the proton 
deformed shell closure. On the other hand, we have mentioned above that 
the deformation energy curve of 32 Mg [see Fig. 3] displays an inflexion point 
at a prolate deformation. This is the first indication of the deformed shell 
effect which determines the shape of the deformation energy curves of the 
heavier isotopes. The values of f3 n and f3 p for this inflexion point are reported 
on Fig. 5. As the i-j/2 shell begins to be filled, the magnitude of the spherical 
neutron shell effect decreases and allows the proton deformed shell effect to 
take over in 36_40 Mg. In these nuclei, one recovers a proton deformation 
similar to that found for 24 Mg. The experimental evidence based on the ex- 
istence of a low energy 2 + state with a large B(E2) value [18] in 32 Mg implies 
that our calculation overestimates the magnitude of the spherical N=20 shell 
effect in this nucleus. With two other Skyrme parametrizations, Sly4 and 
SkP [9], we also find that the deformation energy curve for 32 Mg presents a 
spherical minimum. To our knowledge, this is also the case for all mean-field 
calculations whether they include a treatment of pairing correlations[27] or 
not. Apart from 32 Mg, data on deformation are few with generally large er- 
ror bars. Still they are consistent with values of f3 larger than those we have 
calculated. We believe that this systematic effect reflects the approximate 
character of the the prescription (3) that we use to define j3 v . Indeed, it is 
implicitly based on a sharp edge liquid drop picture which neglegts surface 
effects which could be of importance in such light nuclei. Keeping this re- 
striction in mind, we note that for light isotopes, the values of f3 n and f3 p are 
similar, while they differ for the last three heavy isotopes. With mean-field 
effective forces, differing deformations are more easy to obtain when there is 
a large excess of neutrons over protons. Indeed, for nuclei close to the stabil- 
ity line, the large overlap between the neutron and proton densities and the 
strongly attractive neutron-proton interaction have a tendancy to suppress 
the differences of deformations. This is indeed what we find for the light 
magnesiums. In view of our results, it appears worthwhile to experimentally 
investigate whether heavy Mg isotopes would not be the lightest examples of 
nuclei with differing neutron and proton deformations. 

On Fig.6, one sees that the calculated root-mean-square radii stay within 
the error bars of the measured data. The onset of deformation in 36 Mg 
causes a slight irregularity in the calculated curves. An irregularity of the 
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same magnitude does not seem to be present in the experimental data around 
32 Mg. 

The single particle level energies in the canonical Hartree-Fock-Bogoliubov 
basis are shown on Fig. 7 for 40 Mg. Apart from a global shift in energy, a 
similar diagram is valid for all Mg isotopes. In 40 Mg, at the quadrupole mo- 
ment of the deformed ground state, the neutron fermi level is close to the 
continuum. The large density of levels at the fermi surface provides an expla- 
nation for the smaller neutron deformation (compared to proton). Indeed, 
because it enhances the effect of pairing correlations, it induces a population 
of orbitals with a less deformation driving character than the 716.5/2 orbital. 

Results qualitatively similar to ours have been obtained in the RMF cal- 
culation of Ren et al.[7]. Both 20 Mg and 32 Mg have been found spherical 
by these authors and 36_40 Mg very deformed. However, the quadrupole mo- 
ments of deformed isotopes are significantly larger in the RMF calculation. 
This may be related to the BCS approximation with a constant gap used by 
Ren et al. 

The N = 28 neutron gap does show up as an accident neither in the S 2n 
curve nor in the evolution of the radii. On Fig. 7, one sees that the value of 
the N=28 single particle gap is about 3 MeV in 40 Mg. This is much smaller 
than the gap that one finds in 48 Ca (5MeV) with a calculation using the same 
effective force [31]. This quenching of the N = 28 gap is sufficient to allow 
the deformed proton gap at Z — 12 to win over the spherical shell effect and 
to lead to a deformed nucleus. Similar results have been obtained in a recent 
study of Si and S isotopes [8]. 

Fig. 8 shows the potential energy curve of the nucleus 42 Si (N = 28). 
The same figure displays the curves of the two nuclei 34 Ne and 46 Si which 
according to our method are on the drip line. The energy curves of these 
three isotopes display a soft dependence on the quadrupole moment. This 
suggests that dynamical quadrupole collective effects may play an important 
role. One notes also that, as was the case already for 40 Mg, there is no 
indication of a spherical shell effect at N = 28 for 42 Si. Similar results have 
been obtained by Werner et al[8] with the HF+BCS approximation. 
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4 Conclusion 



In this paper, we have reported the first calculation of the deformation prop- 
erties of drip line nuclei using a Hartree-Fock-Bogoliubov approach in which 
pairing correlations in the continuum are correctly treated. 

Our work shows that deformation effects cannot be neglected in a descrip- 
tion of the structure of neutron rich Mg isotopes. In particular, it provides 
a negative answer to the question raised by Tanihata et al [1] as to whether 
all drip line nuclei are spherical. 

Another conclusion of this work, is that the N=28 spherical shell closure 
is significantly weakened in light nuclei close to the drip line. As a conse- 
quence, we find that the deformation properties of heavy Mg isotopes are 
more influenced by the proton deformed shell effect at Z—12 than by the 
neutron spherical shell closure at N=28. Due to pairing correlations, the 
population of neutron deformation driving orbits favours the occurence of 
deformation in the last bound Mg isotopes. 

Finally, we find that with the prescription that we introduce to define 
deformations, the large excess of neutrons with respect to the protons leads 
non zero differences (f3 n — /3 P ) in the heavy Mg isotopes. A confirmation 
of this property would be the existence of low lying quadrupole isovector 
excitation for isotopes close to the drip lines. 

Therefore, the several original properties predicted by our calculation give 
an additional incentive to the pursuit of a more thorough investigation of this 
section of neutron drip line. 
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Figure Captions 

Fig. 1 Calculated and experimental of Mg isotopes for A = 22—34. with 
two Skyrme force parametrizations: (a) SIII and (b) SLy4. In both 
cases, two strengths of the pairing force are tested: Vq = 1000 and 700 
MeVfm 3 . 

Fig. 2 Calculated and experimental S^n from the proton to the neutron drip 
lines. Our results are obtained with the SIII parametrization. Re- 
sults of RMF and FRDM calculations are taken from refs. [7] and [35], 
respectively. 

Fig. 3 Deformation energy curves calculated for the even 20_40 Mg isotopes 
with SIII as a function of the axial quadrupole moment. The origin of 
the energy is taken at the minimum of each curve. 

Fig. 4 Deformation parameter /3 for protons and neutrons in 40 Mg calculated 
as a function of the cut-off radius r in . See eqs. (3) and (4). 

Fig. 5 Deformation parameter j3 (upper panel) and total quadrupole moments 
(lower panel) of the Mg isotopes versus A. Ground state values are con- 
nected by lines. Isolated symbols correspond to the secondary minima 
found for some isotopes. Black symbols correspond to absolute or sec- 
ondary minima, open symbols to shoulder in the deformation energy 
curves of Fig. 3. The (3 deduced from the experimental data have been 
taken from refs. [36] and [37] for 24 < 26 Mg and from [18] for 32 Mg. 

Fig. 6 Calculated and experimental rms radii versus A. The theoretical values 
are joined by lines. The black stars are the matter rms radii of ref. [24], 
and the white stars are the charge rms radii of ref. [37]. 

Fig. 7 Single-particle energy diagram of 40 Mg as a function of the total quadrupole 
moment. These energies are the diagonal matrix elements of the HF 
Hamiltonian in the canonical basis. The upper and lower parts are for 
neutrons and protons respectively. 

Fig. 8 Deformation energy curves of 34 Ne and 42 > 46 Si as a function of the axial 
quadrupole moment. 
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